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Os| ■ Abstract 
in 
o 

\ We study the effect of random matter density fluctuations in the sun on 

resonant neutrino conversion in matter. We assume no specific mechanism 
O for generation of the fluctuation and we keep the amplitude and correlation 

length as independent parameters. We do not work under the approximation 
that fluctuations have spatial correlations only over distances small compared 
to the neutrino oscillation lengths. Instead we solve numerically the evolution 
equation for the neutrino system including the full effect of the random mat- 
^ ' ter density fluctuations of given amplitude and correlation length. In order 

to establish the possible effect on the MSW solutions to the solar neutrino 
problem we perform a global analysis of all the existing observables including 
the measured total rates as well as the Super-Kamiokande measurement on 
the time dependence of the event rates during the day and night and the 
recoil electron energy spectrum. We find the effects of random noise to be 
larger for small mixing angles and they are mostly important for correlation 
lengths in the range few 100 km < Lq < few 1000 km. They can be under- 
stood as due to a parametric resonance occuring when the phase acquired by 
the oscillating neutrino state on one fluctuation length Lq is a multiple of 2n. 
We find that this resonant parametric condition is mainly achieved for low 
energy neutrinos such as the pp-neutrinos and therefore its effect is mostly 
seen on the total event rates while the other Super-Kamiokande observables 
are very marginally sensitive to the presence of noise due to the higher energy 
threshold. 
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I. INTRODUCTION 



Solar neutrinos were first detected already three decades ago in the Homestake experi- 
ment and from the very beginning it was pointed out the puzzling issue of the deficit in 
the observed rate as compared to the theoretical expectation based on the standard solar 
model with the implicit assumption that neutrinos created in the solar interior reach 
the Earth unchanged, i.e. they are massless and have only standard properties and inter- 
actions. This discrepancy led to a change in the original goal of using solar neutrinos to 
probe the properties of the solar interior towards the study of the properties of the neutrino 
itself and it triggered an intense activity both theoretical as well as experimental, with new 
measurements being proposed in order to address the origin of the deficit. 

On the theoretical side, enormous progress has been done in the improvement of solar 
modelling and calculation of nuclear cross sections. For example, helioseismological ob- 
servations have now established that diffusion is occurring and by now most solar models 
incorporate the effects of helium and heavy element diffusion @,|]]. From the experimental 
point of view the situation is now much richer. Four additional experiments to the original 
Chlorine experiment at Homestake || have also detected solar neutrinos: the radiochemical 
Gallium experiments on pp neutrinos, GALLEX || and SAGE and the water Cerenkov 
detectors Kamiokande [[J and Super-Kamiokande PJTUf. The latter have been able not 
only to confirm the original detection of solar neutrinos at lower rates than predicted by 
standard solar models, but also to demonstrate directly that the neutrinos come from the 
Sun by showing that recoil electrons are scattered in the direction along the Sun-Earth axis. 
Moreover, they have also provided us with useful information on the time dependence of 
the event rates during the day and night, as well as a measurement of the recoil electron 
energy spectrum. After 825 days of operation, Super-Kamiokande has also presented pre- 
liminary results on the seasonal variation of the neutrino event rates, an issue which will 
become important in discriminating the MSW scenario from the possibility of neutrino os- 
cillations in vacuum flTT| , |T2"|| . At the present stage, the quality of the experiments themselves 



and the robustness of the theory give us confidence that in order to describe the data one 
must depart from the Standard Model (SM) of particle physics interactions by endowing 
neutrinos with new properties. In theories beyond the SM, neutrinos may naturally have 
new properties, the most generic of which is the existence of mass. It is undeniable that 
the most popular explanation of the solar neutrino anomaly is in terms of neutrino masses 



and mixing leading to neutrino oscillations either in vacuum |13| or via the matter-enhanced 



MSW mechanism . 

The standard MSW analysis is based on a mean-field treatment of the solar background 
through which the neutrinos propagate. In this approximation the global analysis of the full 



neutrino data sample described above |15j leads to the existence of three allowed regions in 



the Am 2 sin 2 29 parameter space for neutrino oscillations 

• non-adiabatic-matter-enhanced oscillations or small mixing angle (SMA) region with 
Am 2 = (0.4-1) x 10" 5 eV 2 and sin 2 (29) = (1-10) x 10~ 3 , and 

• large mixing (LMA) region Am 2 = (0.2-5) x 10~ 4 eV 2 and sin 2 (26 l ) = 0.6-1. 

• low mass solution (LOW) Am 2 = (0.3-2) x 10~ 7 eV 2 and sin 2 (2fl) = 0.8-1. 
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There are several works in the literature |16|- pl|1 where corrections to such mean-field 
picture have been studied. The influence of periodic matter density fluctuations of given 
amplitude and fixed frequency above the average density on resonant neutrino conversion 



was investigated in Refs. |L6||17|] . In Ref. |T6[ a parametric resonance is found when the fixed 
frequency of the perturbation is close to the neutrino oscillation frequency. This approach 
however gives not answer to the physical origin of such fixed frequency perturbation. 

More recently, the main approach to fluctuations that has been pursued fTffiffiH is to 
model the matter density as a Gaussian random variable (white noise) . In this approach the 
number of free parameters remains the same as in the case of fixed frequency perturbations 
- two, the perturbation amplitude and the correlation length- but the presence of white 
noise in the sun is doubtless since there are many mechanism to generate random density 
perturbations. For technical reasons these analysis where performed under the assumption 
that fluctuations have spatial correlations only over distances small compared to the neutrino 
oscillation lengths. Within this approximation the conclusions obtained were that such 
fluctuations in the solar electron density can significantly modify the MSW solutions to 
the solar neutrino problem (SNP) provided that their relative amplitude near the MSW 
resonance point can be as large as few percent. 

In Ref. [^TJ criticisms to these results were raised based on two facts: i) the unexistence 
of a plausible source for such 5-correlated fluctuations in the vicinity of the MSW resonance 
point and ii) whatever its origin, the effect of the density perturbation is maximum in the 
regime where the short correlation-length approximation fails. In particular in Ref. pl[ 



they concentrate on helioseismological waves as origin of the perturbation and in particular 
on g-waves whose amplitude increases with the solar depth and can, in principle, reach the 
interesting values to affect neutrino propagation. However they conclude that such g-waves 
do not affect the MSW neutrino conversions since the wavelength for the lower modes, 
for which the largest amplitudes are possible, is much longer than the characteristic MSW 
neutrino oscillation length. 



It has been recently argued p2 |, however, that in the magnetohydrodynamical (MHD) 
generalization of the Helioseismology the objection in Ref. does not hold. Assuming 
modest central large-scale magnetic fields (B =1-10-100 Gauss) one can find magneto- 
gravity eigenmodes with much shorter wave lengths for density perturbations Xmhd ~ 
200 — 2000 km |22] that is comparable with the neutrino oscillation length at the MSW 
resonance for large and small mixing angles correspondingly. 

In this paper we revisit the problem of the effect of matter of density fluctuations in the 
sun on the MSW solutions to the SNP. In our approach we assume no specific mechanism 
for generation of the fluctuation and we keep the amplitude and correlation length as inde- 
pendent parameters. There are two main differences in our analysis as compared to those in 
Refs. JL8-21]. First we do not work under the approximation that fluctuations have spatial 
correlations only over distances small compared to the neutrino oscillation lengths. Instead 
we solve numerically the evolution equation for the neutrino system including the full effect 
of the random matter density fluctuations of given amplitude and correlation length. Sec- 
ond, in order to establish the possible effect on the MSW solutions to the SNP we perform 
a global analysis of all the existing observables including not only the measured total rates 
but also the Super-Kamiokande measurement on the time dependence of the event rates 
during the day and night, as well as the recoil electron energy spectrum. In particular we 
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include the regeneration effects when neutrinos cross the Earth [23| which were neglected in 



Ref. [20 



The outline of the paper is as follows. In Sec. [TJ we discuss our approach to the solution 
of the neutrino evolution in the presence of density fluctuations. In Sec. |II A| we briefly 



summarize the standard analytical approach based on the short correlation length approx- 
imation and in Sec. [II B| we discuss our numerical treatment and present our results for 



the survival probabilities as a function of the relevant oscillation and noise parameters. In 
Sec. pTC] we interpret our results for the enhancement of the survival probability in the lan- 



guage of parametric resonance in MSW conversions. Section [TT1] is devoted to the statistical 
analysis of the solar neutrino observables in the framework of the MSW solutions of the SNP 
in presence of the noisy density fluctuations. We study the variation of the allowed regions 
of the SNP for different combinations of observables when noise fluctuations with different 
correlation lengths are included. Our results are summarized in Figs. and Table |I[ We 
show that even for noise levels as large as 4% the relative quality of the three allowed regions 
of the MSW solutions to the SNP depends on the value of the correlation length studied in 
the wide range Lq = 70-10 4 Km and that the three allowed regions of the MSW solutions of 



the SNP remain valid at the present level of solar neutrino experiments. Finally in Sec. [TV 
we discuss our results and summarize our conclusions. 



II. MSW SOLUTIONS TO SNP FOR NOISY MATTER DENSITY IN THE SUN. 

Some mechanisms to rise density perturbations in the Sun were proposed in literature in 
connection with the MSW solution to the SNP. One of them with g-waves in Helioseismology 
as a plausible source of matter noise fails being applied to the MSW neutrino conversions 
since the wavelength of such modes happens to be too long, X g ~ O.li?© ^> l osc to influence 
neutrino oscillations [plj| . This X g is the wave length for low radial degree n < 3 for which 



largest g-mode amplitudes 5p(r)/ p Q (r) ~ 4 % are possible |pl 



However, in the MHD generalization of the Helioseismology such objection [[H]] does not 
apply. Assuming modest central large-scale magnetic fields (B = 1 — 10 — 100 Gauss) one can 
find magneto-gravity eigenmodes with much shorter wave lengths for density perturbations 
Xmhd ~ lose ~ 200 — 2000 km [^2| comparable with the neutrino oscillation length at the 



MSW resonance for large and small mixing angles correspondingly. 

Note that standard Helioseismology corrections to the standard solar model (SSM) (ne- 
glecting magnetic fields in the Sun) give density fluctuations deep in solar interior at a low 
level Sp(r)/po(r) ^ 1% [|4[]. This analysis is done using the solution of the inverted prob- 



lem in Helioseismology based on the corresponding integral equation for which the kernels 



are built on the full set of p-mode waves calculated and observed on the photosphere p5 
However, such analysis does not touch both g-modes that are still invisible on the surface 
of the Sun (from SOHO satellite) and has no relation to MHD modes found in [22 . 



A. Generalized Parke formula for averaged evolution equation 

In this subsection we show that the Schrodinger equation approach for noisy matter EC 



generalized here for an arbitrary matter density perturbation correlator (not only (^-correlator 
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like in POD is equivalent to the Redfield evolution equation leading to the generalized Parke 
formula for the survival probability in the presence of noise ||21|| . 

Let us assume the presence of regular density waves Re fy\? n £p n (£)j exited somehow 

n 

within the solar interior. They could be, for instance, the MHD matter density waves, Sp n (z), 
which appear in the 1-dimensional solar model with the exponential matter background 
profile po(z) ~ exp(— z/H), H ~ O.li? , in the presence of gravity g = (0,0, -g) and an 
external constant magnetic field B = (B ,0,0). Such waves obey the dispersion relation 
u n = u(n, B , k x , k y ) with the periods T = 2ir/uj n ~ few days and they are quite different 
from the g-modes in Helioseismology. In particular, unlike the helioseismological g-modes, 
they have very short wavelength along the z-axis A 2 ~ R & /n <C X g for large node numbers 
n ^> 1 acceptable in the model . 



Assuming these density perturbations added to the SSM background density profile po(t), 
we write the master Schrodinger equation for MSW conversions of two neutrino flavors, 



Vy, 



u e \ = f H e s 2 S 
v y J 1 s 2 5 




(1) 



where in the diagonal entry H e = V ey (t)[l + Re (^2C n (j) n (t)^j] — 2c 2 5 any density eigenmode 

n 

£, n (t) = C n (j) n {t) = 8p n {t)/ po(t) has the small amplitude, | £„(£) |<C 1. c 2 = cos 29, s 2 = sin26 l 
and S = Am 2 /4:E are the neutrino mixing parameters; V ey (t) = G p^/2(p(t) /m p )(l — Y n ) and 
Vey — Ks(^) = G F\^2(p(t) /m p )(l — 3Y n /2) are the neutrino vector potentials in the Sun for 
active-active (y= x) and for active-sterile neutrino conversions correspondingly. They are 
given by the neutron abundance Y n = m p N n (t) / p(t) where for neutral matter the relation 
Y e = Y p = 1 — Y n has been used and the SSM density profile pit) is given by BP98 model. 
In what follows we will discuss only conversion into active neutrinos. 

Using the survival probability P ee = v*v e and the auxiliary functions I = Im (v^y), R = 
Re {y* e Vy) one can derive from the master equation above the equivalent system of dynamical 
equations. After averaging of those dynamical equations over small density perturbations 
I £,n(t) |<S 1, such system takes the form: 




-2n{t) -{V ey {t) - 2Sc 2 ) 
V ey (t) - 2Sc 2 -2n{t) -Ss 2 
+5s 2 



(2) 

with V =< Vii >, X =< X >, 1Z =< 1Z >. This system of equations is similar to Eq. (3.14) 



f n \ 
x 

\V-l/2) 



in 



2U| but here the matter perturbation parameter n(t) is of the form: 

«(*) = 111 I V e y(t)V ey (t 2 )(C n C m )<j) n (t)(j) m (h)dt 2 , (3) 

m,n *l 

where in averaging (...) we use that: (i) C n are uncorrelated random variables which are 
Gaussian distributed with vanishing mean: (C n ) = 0, (ii) different modes (with different 
node number) are uncorrelated, (C n C m ) = A n S nm . 

It is reasonable to assume that even for a regular eigenmode "n" excited somehow within 
the solar interior, different neutrinos emitted from different starting points to in the core 
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propagate to the Earth along parallel rays crossing different profile realizations of the same 
mode £„,(£). The phase of the wave entering in C n is random, (C n ) = 0, as well as a starting 
point to is random for given t). In other words, the averaging over random phases is equiv- 
alent to the averaging over production point distribution (see the discussion in Sec. [II B|) . 

Thus, one can consider the sum of multimode regular perturbations as random density 
perturbations J2C n (fin(t) — > £(£), (£(£)) = 0, for which the 5-correlated noise 

<f(ti)£(* 2 )> = L (^t 1 ) 2 )5(t 1 -t 2 ) 

is the particular case leading from Eq. (||) to n(t) = V 2 y (t)(^ 2 )L /4 f20j . 

We can now identify the entries in the Hamiltonian Eq. (0) with the corresponding 



terms in Eq. (77) derived in Ref. fl9j through the averaging of the Redfield equation for the 
density matrix. One finds full coincidence with the notation in Ref. ||19|| : p\ = X, p 2 = 1Z, 
p 3 = V - 1/2 and a(t) = n{t), -2(M 3 + b) = V ey {t) - 5c 2} M 2 = 0, M x = 5s 2 . 

For slowly varying variables V ey (t) and n(t), one can diagonalize the Hamiltonian in 
Eq. (0) obtaining the generalized Parke formula (Eq. (86) in Ref. [|19j) 

Pee(t) = - + Q - P?j exp (-2 £ jo(x)dx^) cos 2e m (t ) cos 29 m (t) (4) 



where the effect of the density fluctuations is contained in the 70 factor 

4k(x)5 2 sI 
A5 2 s 2 2 + (V ey (x) - 5c 2 y 

Pj is the probability for level crossing as one passes the resonant point PB[ and 6 m is the 
mixing angle in matter. 



*(*) = ,„., 2 v ; 2 ^ ■ (5) 



B. Computer simulation of the master equation 

There are some remarks to the analytic approach shown above that forced us to perform 
a numerical calculation. First, the generalized Parke formula in Eq. ([5]) does not describe 
neutrino conversions for large correlation lengths Lq ^> l OS c in an appropriate way giving a 



huge discrepancy with the results from numerical calculations ^IJ . Alternatively, in Ref. |2T 
the survival probability was evaluated using the standard Parke formula for each neutrino 
ray and then averaging this result over 200 random density profiles of the Cell type of length 
Lq. Below we refer to this procedure as the "Cell" model. Their result for large L , however, 
presents a strong dependence on the resonance position within a Cell. 

In our approach we directly solve the Schroedinger equation for the two-neutrino system 
with the following procedure. Not appealing to any origin of the matter density perturba- 
tions we consider different levels of the random density fluctuations (£(r)) = (Sp(r) / Po( r )) — 
added to the background matter density po( r ) m SSM which we take to be the BP98 density 

p(r) = p«,(r)[l+£(r)] , (6) 



where the parameter £ = y (£,(r) 2 ) measures the amplitude of the perturbation. For a given 
value of £ and correlation length L we generate a density profile of the type in Eq.(0). 
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The function £ (r) is constructed as steps of constant values £_i for each step i of length L 
from the production point to the edge of the Sun. The £j numbers are randomly generated 
following a Gaussian distribution of mean (£(r))=0 and dispersion £. For illustration in 
Fig. [l| we show the numerical profile generated using this procedure for L = 700 km and 
£ = 0.1. 

Substituting this realization of the matter density into the Schrodinger equation for two 
neutrino flavors we have solved the Cauchy problem for different starting points ro within 
the solar core r < R core = 0.3R Q with the same initial condition u e (ro) = 1. The production 
points r are chosen to be placed in knots of a 30 x 60-net that covers the cross-section of 
the core hemisphere, r < 0.3R®. In other words, O.Oli? ~ 7000 km is the cell size of a 
k- rectangle chosen within core. 

In this way, for 30 x 60 = 1800 ro-points, we have obtained a set of the complex wave 
functions, u a (r, r ) =| u a (r,r ) | exp(z$(r, r )), a=e,x, from which one can easily get the 
survival probability at the surface of the Sun, P ee (R Q ,r ) = P ee (R e ,r , s 2 , 5, L , £) and the 
survival probability on the day-side of the Earth after propagation in vacuum through the 
solar wind (e.g. for u e — ► oscillations), 

P£»(r ) = P£»(r , si 5, L , £) = P ee (P , r ) + |[1 - 2P ee (P , r )] - 

--S2C 2 [vl(R Q , r o )^(-R , r ) + ^*(P , r )u e (R Q , r Q )] . (7) 

Then, assuming spherical symmetry, we have averaged the survival probability in Eq. (|7J) 
over the production points ro. The averaging means the multiplication by the weight factor 
defined as the local //-source distribution Si(r k), given by SSM BP98 for each r Q k and each 
neutrino flux type i=pp, Be, pep,..., resulting in: 

i 1800 

(P e t^(4 5, L , £)> = ygoo £ P^MS^) . (8) 

Notice that Eq. fl7|) is equivalent to the standard expression for the day-side survival prob- 
ability P^ y (r ) = P^ un (r )P^ arth + Pg un (r )P£ arth =| c 2 ^ e - csip„ | 2 + | s 2 ^ e + csijj^ | 2 
where the complex wave functions ip a = z/ a (P ,r o ) are given at the surface of the Sun and 
during the day p^ e arth = 1 — p^arth _ cog 2^ During the night, solar neutrinos cross the 



Earth before reaching the detector and regeneration of iVs is possible ||23|| . In order to 
take into account this effect we compute the probability p^ e arth by integrating numerically 
the differential equation that describes the evolution of neutrino flavors in the Earth. This 
probability depends on the amount of Earth matter travelled by the neutrino, or in other 
words, in its arrival direction which is usually parametrized in terms of the zenith angle <3>. 
Thus, in general, the survival probability for a neutrino of given source i arriving at a given 
zenith angle <£> is given by 

(2(P e ^(4 5, L , £)> ~ 1) (s 2 - P 2 E e aHh (4, 8, L , £; $)) 

Such probability remains a function of the fundamental neutrino parameters 5 and s 2 ., as 
well as of two noise parameters £, L Q . 
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In principle we should also average over different realizations of the density perturbations 
with the same level of noise £ and correlation length L . We discuss next that our averaging 
procedure over the grid of production points r is equivalent to the numerical average over 
an ensemble of electron densities. 

Notice that for parallel rays which are directed along the z-axis to the Earth at a fixed 
distance z = z Q from the center, the density profile Eq. @ has different density amplitudes 
since z — t\ only for one ray in equiliptics, for other k rays (k ^ 1) the hypotenuse is 
longer, > z. Thus, rays are not equivalent to each other. This means that considering 
parallel rays and for any k-ray substituting the same final distance z = R & into the wave 
functions u a (R Q , r fc) we automatically took into account different density profile realizations 
including a matter noise. Then integrating (summing) over tq in Eq. (||) we have averaged 
simultaneously over different realizations of noise. 

Alternatively we may think of our averaging procedure in the following way. In our 
distribution of starting points in the grid we have several points j with the same value of ro 
but located at different distances yj from the Sun-Earth z axis. All these points encounter 
different realizations of the matter density in their way as they move in their path parallel to 
the z axis since for each of them rj(t) = \J Zj(t) 2 + y 2 is different and so it is the profile p they 
are subject to at each time t. So our averaging over 1800 initial points can be understood 
as an integration over the starting point r times an average over different matter density 
realizations j for each r . 

In Fig.0 we show the averaged survival probabilities < P e d e ay (s|, 5, L , £) > averaged for 
the 8 B production point distribution, as a function of 5 = Am 2 /4E for two values of the 
mixing angle si = 0.006 (SMA) and s 2 , = .78 (LMA) for level noise £ = 4 % and different 
values of the correlation length L . Also shown in the figure is the corresponding probability 
for the noiseless case. As shown in the figure, even for modest noise level £ = 4 % we find 
relatively large effects. This is specially the case for the SMA value. For LMA the effect is 
large only for short correlation lengths. We next discuss the interpretation of this behaviour. 

C. Parametric resonance for MSW conversions in noisy matter 

In Fig. | we show isolines AP ee = [< P^ y (s 2 2 , S, L , £) > - < P^ y (s 2 2 , S) >] (averaged for 
the 8 B production point distribution) in the plane 5, L for noise level £ = 4 %. P^ y (sl, 5) is 
the survival probability in the absence of noise. Figure [|a corresponds to fixed s 2 = 0.0063 
while Fig. [3|.b corresponds to s| = 0.79. 

One can see in Fig. ^.a a wide spectrum of the domain sizes L , for which the noisy 
solution gives large difference for the v e - suppression, AP ee = 0.3 — 0.5, even for a modest 
noise level (=4 % )!. Largest enhancement occurs for values of 5 = (1-3) x 10~ 6 eV 2 /MeV. 
For characteristic values of Am 2 ~ (0.4-1) x 10~ 5 eV 2 in the SMA region the energy values 
for such enhancement in Fig. |3].a corresponds to E — Am 2 (eV 2 ) / (45) MeV ~ 0.33-2 MeV 
in the interesting range for the existing solar neutrino experiments. 

For the larger mixing angle the enhancement is smaller (see Fig. |3].b). The maximum 
effect of about (~ 0.13) appears somewhere near L$ ~ 70 — 100 km. 

For both cases we can interpret this enhancement as being due to parametric effects 
in neutrino oscillations. The parametric resonance implies a synchronization between the 
system eigen-oscillations having MSW frequencies u = 2tt jl m and the parameter variations 
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given by a changing size of density fluctuations Lq, uj = 2n / Lq. Here the neutrino oscillation 
length in medium l m is given by 

lm{r) = = , (10) 

(cos 29 -h/k{r)) 2 + sin 2 29 



where l v = 4nE/Am 2 is the oscillation length in vacuum, /o( r ) — 2ixm p j \/2GFp(r) is the 
refraction length. 

The synchronization condition (parametric resonance condition) states that the phase 
acquired by the oscillating neutrino state on one fluctuation length L should be a multiple 
of 2tt @ 

dr—— = 2-nk, k = 1,2, 3, ... (11) 
o l m [r) 

One can simplify this condition for SM A MSW when the neutrino energy differs considerably 
from the MSW resonance energy, | cos 29 — I u /Iq |^> sin 29. In this approximation substi- 
tuting the mean density p = p we obtain from Eq. ( ]TT| ) the simple formula for parametric 
resonance in the case of SMA MSW oscillations 



. COS26 1 1 . k . , 

l -^~T l= T- (12) 

According to Eq. (0), two types of resonance regimes may appear depending on the 
value of L : 

• A resonance at low energy when 5 5r = GfPr/ {\f2m p cos 29) (p < pr). In this 
case, for Lq <^Iq, one obtains the parametric resonance condition l m ~ l v ~ Lq. 

• A resonance at energy well above the MSW resonance, 5 C ^ (p > pn). In this case 
the parametric resonance condition is l m ~ Iq ~ Lq. 



A more interesting case occurs when the synchronization condition in Eq. (|TTD is achieved 
in the proximity of the MSW resonance (5 = 5r). In this case we can substitute V^f = 
l v j sin26*, and still assuming a short correlation length, Lq < Ar, compared to the thickness 
of the resonant layer Ar, we get much longer Lo-sizes 

Lo = kl - = T^V9- (13) 

Or sin 20 



The parametric resonance is stronger for this case 

Since for SMA MSW conversions the thickness of the resonant layer is Ar = 
p(dp/ dr)~ l / cot 2$ ~ O.l-R / cot 29 ~ 7000 km the condition in Eq. (|13|) remains valid for 
some first numbers k = 1,2,. ... In other words, the condition 7000 km > Lq = kir/Sjism29 
implies a lower bound on the possible values of 5r for which the parametric resonance is 
possible. For instance for the mixing angle in Fig. ^.a sin(2#) = 0.08, the value Sr ~ 10~ 6 
eV 2 /MeV is only marginally allowed for k = 1. 

Conversely from the MSW condition V£ s = 250 km(E/MeV)/(Am 2 /10^ 5 eV 2 )/sin2#) 
substituting the SMA MSW parameter values Am 2 ~ 10~ 5 eV 2 , sin 29 ~ 0.1 we find 
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1™ S ' SMA ~ 2500(i?/MeV) km. Moreover we see in Fig. ^.a that domain sizes Lq = 
700 — 2000 km are appropriate for maximum enhancement which implies that the MSW 
resonant parametric condition Eq. (O) is mainly achieved for k = 1 in the case of low 
neutrino energies like the pp-neutrinos seen in GALLEX and SAGE. 

Let us note that the parametric resonant condition in Eq. (O) is rather opposite to what 



III ; 



was assumed in Eq. (3.18) p0[ , L = 0.1 l m . In other words, 5-correlated noise, L <C I 
is an analytic approximation which is appropriate only for small values of the correlation 
length. To illustrate this we plot in Fig. [| the survival probability (P ee (s§, 6, Lq, £)) as the 
function of L for different values of the level noise and fixed 5 = 10~ 6 eV 2 /MeV. Figure |].a 
corresponds to s 2 2 = 0.01 (1™ S ~ 7000 km) and Fig. |.b to s\ = 0.7 (1™ S ~ 700 km). 

One can see in Fig. f| different heights of the "bumps" of the survival probability for 
different noise levels whereas the position of the peaks are similar cutting sharply at the 
right edge L = l^ s - position of the strongest parametric resonance Eq. ([H|) for k — 1. 
If 8 7^ 5r, or energy E is far from the MSW resonance value the parametric resonance 
still takes place but from corresponding Eq.fljjD we find that lower L ~ l m < 1™ S are 
appropriate. Only if L <C l m the 5-correlated regime starts. 

It is interesting to compare our results in Fig. ^.a with the numerical results obtained for 
the "Cell" model plotted in Fig. 1 of Ref. pTj] where authors averaged the ordinary Parke 
formula over random ensemble of 200 density profiles of Cell type (dashed, dot-dashed and 
solid thick lines on that figure). These 200 profiles correspond to 200 neutrino creation sites 
within the core only in contrast with 1800 points in our direct numerical simulation of the 
Schrddinger equation Eq. (JJ). 

Notice that in that figure 5 = 26 Fig. |- We find that both figures present similar behaviour 
in the short correlation length regime and a comparable enhancement for L Q ~ 10 2 -10 4 
which ends at L ~ 1™ S . But for L Q 3> 1™ S the results in Fig. 1 of Ref. |2T| present a strong 
dependence on the position of the resonance within the cell. The result of our numerical 
calculation is closer to the "Cell" model result with randomly distributed cell positions 
(dot-dashed line in Fig. 1 of Ref. f2~Ifl ) and the mean (P ee (L )) tends to the noiseless MSW 
survival probability as expected for large L . 

Concluding this section we want to remark that the numerical approach presented here 
with the averaging of the solutions P ee (r ) over noise realizations is quite different from any 
previous approaches with the averaging of the Schrddinger equation itself before obtaining 
of a solution. Only the straightforward numerical solution of the Schrddinger equation is an 
appropriate way to tackle the problem in the general case of arbitrary correlation lengths. 



III. FITS: RESULTS 

A. Data and Techniques 

In order to study the possible values of neutrino masses and mixing for the oscillation 
solution of the solar neutrino problem in the presence of noisy matter density fluctuations in 
the sun, we have used data on the total event rates measured in the Chlorine experiment at 
Homestake ||, in the two Gallium experiments GALLEX and SAGE ||f7| and in the water 
Cerenkov detectors Kamiokande and Super-Kamiokande shown in Table |5|. Apart from the 
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total event rates, we have in this last case the zenith angle distribution of the events and the 
electron recoil energy spectrum, all measured with their recent 825-day data sample [[nj. 

For the calculation of the theoretical expectations we use the BP98 standard solar model 
of Ref . [ 29 1 . The general expression of the expected event rate in the presence of oscillations 
in experiment i is given by Rf 1 : 



Rf= E <f>h fdE u X k (E u ) x [a e>l (E u )(P e k e (s 2 2 ,5,L ,O) 

fc=l,8 



(14) 



+a x , i (E v )[l-{P? e (sl5,L ,£)) 



where E v is the neutrino energy, is the total neutrino flux and is the neutrino energy 
spectrum (normalised to 1) from the solar nuclear reaction k with the normalization given 
in Ref. p9[ . Here a e ^ (cr x ,i) is the u e (i> x , x = /i, r) interaction cross section in the Stan- 
dard Model with the target corresponding to experiment i. For the Chlorine and Gallium 
experiments we use improved cross sections a e ^(E) from Ref. [30|. For the Kamiokande and 



Super-Kamiokande experiment we calculate the expected signal with the corrected cross 
section as explained below. (-P e fc e ) is is the yearly averaged v e survival probability at the 
detector, as given by Eq. (||) after averaging over arrival directions $. Note that (-P e fc e ) is a 
function of the oscillation parameters as well as the noise parameters. 

We have also included in the fit the experimental results from the Super-Kamiokande 
Collaboration on the zenith angle distribution of events taken on 5 night periods and the day 
averaged value |10| . We compute the expected event rate in the period a in the presence of 



MSW oscillations as, 



R th 



,,,, - a— / , dr £ fc \dE v \ k (E u ) x a e>sk (E v )(P^(sl5,L ^;r)) (15) 

^■T a J T (COS ^ minia ) k = l 8 

k I „2 



+<T XtBk (E v ) l-(P e fc e (45,L ,e;r)) 



where r measures the yearly averaged length of the period a normalized to 1, so Ar a = 
r(cos $ m ax,a) — t(cos $> m in,a) = -500, .086, .091, .113, .111, .099 for the day and five night 
periods. Notice that the dependence of (P^s^S, L ,£;t)) on r comes only from the de- 
pendence of the Earth regeneration probability p^ e arth on the different Earth matter profile 
crossed by the neutrino during the five night periods. 

The Super-Kamiokande Collaboration has also measured the recoil electron energy spec- 
trum. In their published analysis [[J after 504 days of operation they present their results 
for energies above 6.5 MeV using the Low Energy (LE) analysis in which the recoil en- 
ergy spectrum is divided into 16 bins, 15 bins of 0.5 MeV energy width and the last bin 
containing all events with energy in the range 14 MeV to 20 MeV. Below 6.5 MeV the back- 
ground of the LE analysis increases very fast as the energy decreases. Super-Kamiokande 
has designed a new Super Low Energy (SLE) analysis in order to reject this background 
more efficiently so as to be able to lower their threshold down to 5.5 MeV. In their 825-day 
data 
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they have used the SLE method and they present results for two additional bins 
with energies between 5.5 MeV and 6.5 MeV. In our study we use the experimental results 
from the Super-Kamiokande Collaboration on the recoil electron spectrum divided in 18 
energy bins, including the results from the LE analysis for the 16 bins above 6.5 MeV and 
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the results from the SLE analysis for the two low energy bins below 6.5 MeV. The general 
expression of the expected rate in a bin in the presence of oscillations, R , is similar to that 
in Eq. (|I~4"D, with the substitution of the cross sections with the corresponding differential 
cross sections folded with the finite energy resolution function of the detector and integrated 
over the electron recoil energy interval of the bin, T min <T< T max : 



a, 



a,sk 



UT \ TAwZ dT'Res{T, T') da «^ gO 
Jo dT' 



The resolution function Res(T, T") is of the form |9|,j3l|1 



Res(T, T') 



27r(0.47,/T / (MeV)) 



exp 



(T-T 



l\2 



0.44T'(MeV) 



(16) 



(17) 



and we take the differential cross section da a {E vl T')/dT' from ]32[ 



In the statistical treatment of all these data we perform a x analysis for the different 



sets of data, following closely the analysis of Ref. with the updated uncertainties given 
in Refs. |28H30f , as discussed in Ref. Hy|. We thus define a \ 2 function for the three set of 
observables x? ates , xLith, and X 2 pectrum where in both xLith, and ^ pectrum w e allow for a 
free normalization in order to avoid double-counting with the data on the total event rate 
which is already included in Xrates- m the combinations of observables we define the \ 2 °f 
the combination as the sum of the different x 2 ' s - 111 principle such analysis should be taken 
with a grain of salt as these pieces of information are not fully independent; in fact, they 
are just different projections of the double differential spectrum of events as a function of 
time and energy. Thus, in our combination we are neglecting possible correlations between 
the uncertainties in the energy and time dependence of the event rates. 



B. MSW Regions 

We present next the results of the allowed regions in the two-parameter space Am 2 , 
sin 2 {26) for the analysis of the different combination of observables. In building these regions, 
for a given set of observables and a certain value of the noise parameters L and £ (in what 
follows we present the results for £ = 4% which is a reasonable large value for the noise 
amplitude) we compute for any point in the parameter space of two-neutrino oscillations 
the expected values of the observables and with those and the corresponding uncertainties 
we construct the function x 2 (Am 2 , sin 2 (2#); Lq, £) bs- We find its minimum in the full two- 
dimensional space of MSW oscillations. The allowed regions for a given CL are then defined 
as the set of points satisfying the condition: 

X 2 (Am 2 ,sin 2 (2^);L ,Oo bs -xLn, of , s (^o,0 < A X 2 (CL, 2 dof) (18) 

where, for instance A% 2 (CL, 2 dof)=4.6, 6.1, and 9.2 for CL=90, 95, and 99 % respectively. 
In Figs. |^-|8| we plot the corresponding allowed regions for £ = 4 % and for five values of the 
correlation length Lq =70, 200, 700, 2000, and 10 4 km. For comparison in the last panel we 
also show the regions in the absence of random noise. 

Figure [5] shows the results of the fit to the observed total rates only. We see in the 
figure that for any value of the correlation length we always find the three allowed regions, 
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SMA, LMA and LOW as in the standard noiseless MSW analysis although their extend in 
Am 2 and sin 2 (29) varies with the value of the correlation length L . The presence of noise 
modifies the shape and size of the allowed regions through two different although related 
effects. Due to the modification of the shape of the electron survival probability, the value 
of the expected rates for a given point on the MSW plane are different once the noise is 
included what modifies the value of x 2 f° r that given point. This leaves also to a shift on 
the value of Xmin use d to define the regions. 

In Table || we show the values of the local Xmin i n the three regions for the different 
values of the correlation length. First thing we notice is that for the analysis of the rates 
only in the LMA region the value of Xmin * s always lower in the presence of noise. This 
can be understood by looking at Fig. p|.b where we see that noise increases the survival 
probability for 2 x 10 -7 ^ 5 ^ 3 x 10~ 6 , ie for lower values of the neutrino energies relevant 
for the Gallium experiments. This increases the rate at the Gallium experiments which is 
underestimated in the standard LMA solution. For the same reason the LMA become larger 
and shifted to lower Am 2 values and it extends to smaller angles. For larger L , noiseless 
LMA solution is obtained. 

For the SMA solution we find that unless for very short or very long correlation lengths, 
for which Xmin * s almost the same as in the absence of noise, for 100 km < L < few 1000 
km, the presence of noise leads to an increase on the value of Xmim or m °t ner words to 
a worse description of the data in the SMA region. In Fig. |awe see that the presence 
of noise leads to an increase on the survival probability for 2 x 10~ 7 ^$ 5 ^ 3 x 10 -6 , so 
there is not enough suppression of the Berilium neutrinos in the relevant Am 2 values of 
the SMA solution. For this reason the rate at Chlorine experiment is increased and the fit 
worsens as compared to the noiseless case. The worsening is maximal for L of few thousand 
km when the effect of noise is maximally resonantly enhanced as discussed in Sec. [11 C . In 
this case the LMA gives a better description of the measured rates. Also the SMA region is 
shifted towards larger values of the mixing angle. Finally we notice that the LOW solution is 
basically unmodified by the presence of noise as it mainly corresponds to values of 5 ^$ 10~ 7 
which are little affected by the presence of perturbations as seen in Fig. 0.b. 



C. Zenith angle Dependence and Super Kamiokande Spectrum 

Figure |6| shows the regions allowed by the fit of both total rates and the Super- 
Kamiokande zenith angular distribution. Also plotted is the excluded region at 99 % CL 
from the zenith angular measurement. As seen in the figure the shape of the excluded region 
is very little dependent of L , as expected. The zenith angular dependence measures the 
regeneration effect on the neutrino survival probability when the neutrinos travel through 
the Earth matter and it is expected to be independent on the details of the sun matter 
modelling. The main effect of the inclusion of the day-night variation data is to cut down 
the lower part of the LMA region for any value of L . Since for short correlation lengths 
the LMA region had been shifted towards lower Am 2 values, the inclusion of the zenith 
angle distribution data leads to a reduction of the size of the LMA region for short Lq. In 
this way, for instance, for Lq = 70 km the LMA region at 99 % CL extends only in the 
range 1.5 < (Am 2 /10~ 5 ) eV 2 < 7- to be compared with 1.5 < (Am 2 /10~ 5 ) eV 2 < 100 
in the absence of noise-. The SMA is also reduced in size as compared to the noiseless 
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case. We also find that once the zenith angle information is included the LMA becomes a 
better solution (lower xLin) f° r -^0=700-2000 km as seen in Table f|. For these intermediate 
L ~ 700 — 2000 km, when the LMA region obtained from the analysis of the rates extended 
to smaller angles, some tiny regions in between the LMA and the SMA are still allowed at 
the 99 % CL. 

In Fig. |6] we show the regions allowed by the fit of both total rates and the Super- 
Kamiokande recoil electron energy spectrum. Also plotted is the excluded region at 99 % 
CL from the spectrum measurement. The main effect of the inclusion of the spectrum data 
is to improve the quality of the LMA solution as compared to the SMA. Comparing with 
Fig. |5| we observe that the regions become larger after the inclusion of the spectrum data. In 
particular the LMA region extends to larger values of Am 2 . This behaviour is also observed 
in the absence of noise. This is mainly due to the flattening of the x 2 function after the 
inclusion of the spectrum data and it is independent on the presence of noise. Notice that the 
larger sensitivity of the mean (P ee ) to the presence of noise occurs for low energy neutrinos 
as discussed in Sec. |J. Having a threshold energy, T th = 5.5 MeV, Super-Kamiokande is 
insensitive to such low energy effects. 

D. Global 

Figure |8| displays the results for the allowed regions from the global analysis of the solar 
neutrino data including the data on the total event rates, the zenith angular dependence 
and the recoil electron energy spectrum. We find that, after all the observables are included 
the three regions remain valid. The shape and size of final LMA solution is very little 
affected by the presence of noise. The SMA solution is maximally deformed for correlation 
lengths few 100 km < Lq < few 1000 km. It is for these values also that the SMA gives 
a worse description of the observables. This behaviour is mainly driven by the effect on 
the total event rates which is the observable more sensitive to the presence of noise. As 
discussed above both the zenith angle distribution data and the spectrum data have very 
little discriminating power on the noise parameters. One also finds that the tiny 99% CL 
"islands" in between the SMA and LMA are allowed for L ~ 700 - 2000 km. 

Finally we notice that the LOW solution is basically unmodified by the presence of 
the noise. For the neutrino energies accessible at existing experiments, the LOW region 
corresponds to values of 5 ^ 10~ 7 which are little affected by the presence of noise. 

IV. DISCUSSION 

In this paper we have studied the effect of density matter fluctuations in the sun on 
the MSW solutions to the SNP. Assuming no specific mechanism for generation of the 
fluctuations we have kept the amplitude and correlation length as independent parameters. 

Our analysis is performed under no assumption on the relative size of the correlation 
length of fluctuations as compared to the neutrino oscillation length. To perform such a 
study we have solved numerically the evolution equation for the neutrino system including 
the full effect of the random matter density fluctuations of given amplitude and correlation 
length. Our procedure is to generate a realization of the density profile for given values of the 
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perturbation amplitude and correlation length and then to solve numerically the evolution 
equation for the neutrino states for that given realization of the density profile and different 
neutrino production points and finally to average the obtained survival probability over 
different density realizations (with the same amplitude and correlation length) and neutrino 
production points. This numerical approach of averaging the solutions P e e(?"o) over noise 
realizations is different from any previous approaches with the averaging of the Schrodinger 
equation itself before obtaining the solution. 

Our results for the survival probabilities are presented in Figs. |2|-f|. We find that the 
effects are larger for small mixing angles. The larger the mixing angle the shorter the 
correlation length needed to observe an effect. For the SMA the larger effects occur for 
correlation lengths in the range few 100 km < Lq ^ few 1000 km. They can be understood as 
due to a parametric resonance occuring when the phase acquired by the oscillating neutrino 
state on one fluctuation length L is a multiple of 2tt. This resonance is maximal when this 
condition is verified close to the MSW resonance. We find that this resonant parametric 
condition is mainly achieved for low neutrino energies such as the pp-neutrinos seen in 
GALLEX and SAGE. 

Next, in order to establish the possible effect of the presence of noise on the MSW 
solutions to the SNP we have performed a global analysis of all the existing observables 
including not only the measured total rates but also the Super-Kamiokande measurement 
on the time dependence of the event rates during the day and night, as well as the recoil 
electron energy spectrum. The result of such analysis is presented in Figs. [5]^ where we plot 
the allowed regions for MSW neutrino oscillations in the framework of two-neutrino mixing 
with the Sun density profile generated from the BP98, after including random noise with 
amplitude f = 4% and different correlation lengths L (70, 200, 700, 2000 y 10000 km). The 
main conclusions are that the total rates are the most sensitive observables to the presence 
of noise. On the other hand when the many degrees of freedom corresponding to the Super- 
Kamiokande spectrum are included the dependence of the allowed mixing parameters on the 
matter noise is smoothed. This is caused by the larger sensitivity of the mean (P ee ) to the 
noise for low energy neutrinos. Due to its higher energy threshold, the Super-Kamiokande 
experiment is mostly insensitive to these effects. For the same reason one expects that the 
Borexino experiment would be more suitable to place bounds both on the level of neutrino 
noise a/< £ > 2 and on the correlation length L . 
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TABLES 



Experiment 


Rate 


Ref. 


Units 


i? BP98 


Homestake 
GALLEX + SAGE 
Kamiokande 
Sup er-K amiokande 


2.56 ±0.23 
72.3 ±5.6 
2.80 ± 0.38 
2.45 ± 0.08 


§ 

1 

[10| 


SNU 
SNU 

10 6 cm~2 s" 1 
10 6 cm" 2 s" 1 


7.8 ± 1.1 
130 ±7 
5.2 ±0.9 
5.2 ±0.9 



TABLE I. Measured rates for the Chlorine, Gallium, Kamiokande and Super-Kamiokande 
experiments. 



Length 


V 2 




x 2 ■ 


i,LMA 






x 2 ■ 


,LOW 






Rates 


Rates 


Rates 


Global 


Rates 


Rates 


Rates 


Global 


Rates 


Rates 


Rates 


Global 






±Zen 


±Spec 






±Zen 


±Spec 






±Zen 


±Spec 




noiseless 


0.4 


6.2 


22.3 


29.2 


2.9 


7.5 


22.1 


27.9 


7.4 


12.5 


26.9 


32.0 


L = 70 km 


0.2 


4.6 


23.4 


28.4 


1.9 


7.1 


21.9 


28.7 


7.4 


12.5 


26.9 


32.0 


L = 200 km 


0.6 


5.5 


23.7 


28.9 


1.5 


7.0 


21.6 


27.6 


7.4 


12.5 


26.9 


32.0 


L = 700 km 


1.5 


8.0 


22.6 


29.8 


1.7 


7.4 


21.0 


27.6 


7.4 


12.5 


26.9 


32.0 


L = 2000 km 


2.8 


7.5 


25.4 


30.5 


2.1 


7.4 


21.6 


27.6 


7.4 


12.5 


26.9 


32.0 


L = 10000 km 


0.2 


4.8 


22.3 


27.6 


2.5 


7.5 


22.2 


27.8 


7.4 


12.5 


26.9 


32.0 



TABLE II. xlain for SMA > LMA and L0W solutions for the total event rates. 
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FIGURES 



150 




FIG. 1. Numerical density profile as a fuction of the radial distance for Lq = 700 km and 
noiselevel y/ (£) 2 = 0.1. 
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(5=Am 2 /4E (eV 2 /MeV) 



FIG. 2. Electron neutrino survival probability for 8 B neutrinos as a function of S = Am 2 /(4.E) 
for two values of the mixing angle in the SMA region (a) sin 2 (26* ) = 0.006 and in the LMA region 
(b) sin 2 (26*) = 0.78 and for several values of the correlation length as label in the figure. In both 
panels the level of noise £ = 4%. 
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FIG. 3. Isocontours of constant [< P^ v {Lq,5) > — -P e d e ay (<5)] values for 8 B neutrinos in the 
ane 5, Lq for noise level £ = 4 % and for two values of the mixing ang le (a) si = 0.0063 and 



si = 0.79. 
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FIG. 4. Electron neutrino survival probability for 8 B neutrinos as a function of the correlation 
length Lq for 6 = 3.4 x 10~ 6 for two values of the mixing angle (a) sin 2 (29) = 0.01 and (b) 
sin 2 (26») = 0.7 and for level of noise £ =2 %, %, 8 %, 10 %. 
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FIG. 5. Allowed regions in Am 2 and sin 2 (2#) from the measurements of the total event rates at 
Chlorine, Gallium and Super-Kamiokande (825-day data sample) for different correlation lengths 
and for a noise level £ =4 %. For the sake of comparison the results for the standard noisless 
analysis is also shown, he darker (lighter) areas indicate the 99% (90%) CL regions. The best-fit 
point used to defined the regions are indicated by a star. 
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FIG. 6. Same as Fig. |5| but including also the data on the zenith angle distribution observed 
in Super-Kamiokande. The shadowed area represents the region excluded by the zenith angle 
distribution data at 99% CL. 
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FIG. 7. Same as Fig. |^ but including also the data on the recoil electron energy spectrum 
observed in Super-Kamiokande. The shadowed area represents the region excluded by spectrum 
data at 99% CL. 
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FIG. 8. Same as Fig. ^ for a global fit to the solar neutrino data. 
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